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COMPUTER PROGRAM TO PREDICT TOE NEWTONIAN 


AERODYNAMICS OF GENERAL BODIES APPROXIMATED BY FIAT PLATES- 
By Ralph E. Graham, Robert H. Lamb, and Paul 0. Rome re 


SUMMARY 


Equations for computing the Newtonian aerodynamics of a body com- 
posed of flat plate segments have been derived and incorporated into a 
digital computer program. A complete description of the computer pro- 
gram is included, and the use of the program is demonstrated by pre- 
dicting the Newtonian aerodynamics of an example configuration consisting 
of four flat plates. 


INTRODUCTION 


The hypersonic aerodynamic? of a three-dimensional shape can be 
approximated by Newtonian impact theory if tne body is divided into a 
sufficient number of flat plates. The purpose of this report is to 
present the derivation of the equations defining the Newtonian aero- 
dynamics of a general body composed of flat plate segments and to present 
a digital computer program which utilizes these equations. No attempt 
is made to compute the effects of the shading of a segment due to the 
presence of an upstream segment. 
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SYMBOLS 


a, b, c angles between the unit normal and the X-, Y-, 
and Z-axis, respectively 


A ; B, C cosine of angles a, b, and c, respectively 


axial-force coefficient, 


-F, 


qS 


ref 


drag coefficient, D 

qS 


ref 


m 


lift coefficient, L 


qS 


ref 


rolling-moment coefficient, 


**X 


* S ref» 


pitching-moment coefficient. 




^refS 


normal-force coefficient. 


-F. 


qS 


ref 


yawing-moment coefficient, 


* S ref E 


pressure coefficient. 


P - P_ 


side-force coefficient, Y 


qS 


ref 


d, e, f angles between the projection of the unit normal into 

the YZ plane and the Y axis , the XY plane and the 
X axis, and the XZ plane and the X axis respectively 
measured from the axis to the projection 


reference length 


Dx, Dy, Dz cosine of angles a, b, and c, respectively 
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I 


i, J, k 
J, K, L 


L/D 

"x 


drag force 
lift force 
force along X-axis 
force along Y-axis 
force along Z-axis 

unit coordinate vectors 

projection of the unit normal into the XY, XZ, 
and ZY planes* respectively 

lift-drag ratio, C T /C r . 

rolling moment 

pitching moment 

yawing moment 



P 

Poo 

q 

s 


inward directed unit vector normal to surface 

local pressure 
free -stream pressure 

free-stream dynamic pressure 
area 


V 

V 

X, Y, Z 
AX, AY, AZ 

a 


unit wind vector 

Cartesian body coordinate axes 

transfer distances along X-, Y-, and Z- axis, 
respectively 

angle of attack, deg 

angle of sideslip, deg 
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H angle between unit normal and unit wind vectors, deg 

0 angle of roll, deg 

Subscripts: 

loc local area of each flat plate 

ref reference area 


DERIVATION OF EQUATIONS 


In order to determine the Newtonian aerodynamics of a flat plate 
segment, the angle tj between the normal to the segment and the 
relative wind must be established. It has been found that the best way 
to determine this orientation is to use the angles which the normal to 
the segment makes in the principal planes (XY, XZ, and YZ planes). 
Figure 1 presents a unit normal vector at arbitrary angles to the prin- 
cipal axes and principal planes. 

From a three view drawing of a flat plate segment any two of the 
three angles d, e, and f may always be determined. From figure 1, 

since tan d = ^: tan e = ■?: and tan f = — then: 

B' A' A 

C _ C/A 
B B/A 


or 


tan d. 


tan f 
tan e 


Thus, knowing any two angles, the third angle may be determined. 


The vector equation of the unit normal may be written 
7 = Dxi + Dyj + Dzk where Dx, Dy, and Dz are the direction cosines. 



From figure 1, Dx thus represents the cosine of angle a, Dy the cosine 
of b, and Dz the cosine of angle c. It now remains to determine 
angles a, b, and c knowing angles d, e, and f. Again from figure 1 



therefore 


but since 


or 


also since 


— = cos e 
J 

J = sin c 
A = sin c cos e 

B . 

— = tan e 
A 


B 

tan e 
sin c = 


sin c cos e 
cos b 

cos e tan e 


sin e 


cos b 
sin c 


C 

B 


= tan d = 


cos c 
cos b 


sin e 


cos c 

sin c tan d 



tan c 


1 

sin e tan d 


tan f 


cos c 
cos a 


cos a 


cos c 
tan f 


similarly 


, cos c 

cos b = — 

tan d 


Thus all the direction cosines have been determined in terras of 
the angles the projections of the unit normal into the three principal 
planes make with the principal axes. It should be noted that the above 






i.fiM rftv ' 
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relations are valid only for angles l°ss than or equal to 90°. If the 
unit normal lies in other octants, the correct signs must he given to 
the cosine of angles a, b, and c. Figure 2 shows the numbering of 
the octants along with the proper signs for the direction cosines of 
a unit normal vector lying in any octant and directed towards the origin. 

Die vector equation of the relative wind at the Bi’ler angles of 
a, 3, and f to the axes system may be written after inspection of 
figure 3 . 


— - = - cos a cos pi - (sin a cos 8 sin ft + sin 3 cos f ) j 

|V | ^ 

+ (sin 3 sin f - sin a cos 3 co3 f) k 

Die normal force coefficient of a flat plate is defined by New 
tonian theory as 

2 

C„- = 2 COS T] 


Die angle rj may be found by taking the dot product between the 
unit normal vector and the unit velocity vector 


Dius 




cos T) = COS T] 


cos t) = Dx(-ccs a cos 3 ) - Dy(sin a cos 3 sin f + sin 3 cos f) 

+ Dz(sin 3 sin f - sin a cos 3 cos f) 

Die contribution of a segment to the of the total body can 
be expressed as 

C N = 2 erz 2 r, (-Dz) 


where S n is the local area (segment area) and S _ is the reference 
loc ret 

area. (See figure k for coefficient sign convention. ) Although the 

and of a flat plate are zero in Newtonian theory, the flat plate 

Cy may be resolved into the body axis system to produce Doth coefficients. 
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and 


_ , 2 loc / _ v 

C A = 2 cog t] g (-Dx) 

ref 


C y = 2 cos" n (l}y) 

ref 


The moment coefficients may be written: 


_ _ AX „ AZ 

C m " C N D “ C A D 


C 


1 


- -C ^ 
r D 


H 


AY 

D 


C = 
n 


p ^ 

S' D 


+ C. 


AY 

D 


where the signs of the transfer distances (AX, AY, and AZ) are taken 
as positive when the direction from the body moment reference center 
to the segment centroid is in the positive axis direction. The re- 
ference length is D. 


APPLICATION OF THEORY 


In applying the theory, consideration must be given to the -problem 
of shading. In Newtonian theory, it is assumed that the flow travels 
a straight pa oh; there is no curvature of the flew around a body. Thus, 
it is possible for a body segment to partially or completely shield or 
shade a down stream body segment from the flow, resulting in zero force 
on the shielded segment. Since no provision has been made in the pre- 
vious relations to account for this shading, each segment must be con- 
sidered in the presence of the total body for each attitude the body 
makes with the wind to determine if shading occurs. 

In approximating the surface of the desired body by using flat 
plate segments, an important aspect to be considered is the local curva- 
ture of the surface. If the direction of the normal to the given sur- 
face segment to be replaced by a flat plate varies a great deal then the 
aerodynamics can be very incorrect. One example of the former is an 
attempt to replace a spherical segment by pie -shaped flat plates. The 
direction of the normal to the spherical segment pie -shaped wedge varies 
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considerably across the shape. The unit normal for the flat plate wedge 
yields a poor approximation to the actual aerodynamics. The correct 
approach is to choose the flat plate segments small enough so that no 
large change in normal direction will occur. 


EXAMPLE 


An example configuration consisting of four flat plates is pre- 
sented in figure 5 to demonstrate the program's use. Included in the 
figure are the unit normal projections for one of the flat plate seg- 
ments to illustrate the methods of measuring the angles d, e, and f. 
The predicted aerodynamics of the example configuration are presented 
in figure 6. 
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APPENDIX 


Program Description 

The method set forth on the previous pages has been programed for 
a digital computer so that the Newtonian aerodynamics of a flat plate 
or a group of flat plates which approximate a general configuration may 
be computed. The shading or flow - see boundary due to Newtonian theory 
is calculated internally by the program. This is accomplished by test- 
ing the value of cos q for each flat plate. If this quantity is zero 
or negative, the coefficients for this flat plate are set equal to zero 
since the plate is shaded. Shading caused by upstream portions of the 
body must be handled externally to the program; that is, the section of 
the body which would be shaded by a protuberance is not read into the 
program. It should be noted that as either a, 3, or 0 change, this 
shading effect will be different. 

For the convenience of the user, several options or unique features 
are provided and are as follows: 

1 The user may either obtain the aerodynamics for each flat plate 
as well as for the total vehicle or may obtain the aerodynamics of only 
the total vehicle. 

2. The user may vary the angles of attack, sideslip, or roll with- 
out restriction. There is no limit to the number of combinations that 
may be used. 

3. The user may make configuration changes to the configuration 
without having to run a new case. An example of such a feature would 
be a configuration having flaps. The configuration could be run with 
one deflection through a, 3, 0 combinations; then the flap could be 
deflected to a new angle and the resulting configuration could be run 
back through all the same a, 3, 0 combinations. 

4. The program also calculates six component derivatives with 

respect to either a r 3* This is accomplished by perturbing either 
a or 9 1° and assu linearity. 

5. If the vehicle is symmetrical about the XZ plane and is ax 
zero angles of sideslip and roll, only half of the vehicle need be read 
in. The program will calculate values for the total vehicle. 
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Program Input 

The following is a detailed description of the input cards and 
the options that may be used. Presented in table I is an example case 
of input for the configuration presented in figure 5* 

Step 1 Title card 

This is a hollerith card which may contain from 1 to 50 characters 
which describe the configuration being run. 

Step 2 SREF, DX Format 2F10. 0 

SREF = Reference area of configuration 
DX = Reference length of configuration 


Step 3 STPCO Format F10. 0 

STPCO = Pressure coefficient 
Step k J, JS, JL, JA, JB Format 51^ 


J = Total number of flat plate segments 
JS Alievs segment printout if 1 

Allows only total configuration printout if 0 
JL = Number of configuration changes desired 
JA = Number of cards contained in step 6 
JB If 1, it is a symmetrical body and only half of the 
configuration is read in, but the total aerodynamics 
are calculated 

If blank, total body is read in 
NOTE: If JB Is set equal to 1, 3 and $ must equal 0 

Step 5 

E(l) D(l) F(l) X(l) Y(l) Z(l) SLOC(l) KKK(l) 


E(J) D(J) F(J) X(J) Y(J) Z(J) SLOC(j) KKK(j) 

E(l) - Angle measured in the XY plane from the X axis toward the 
Y axis 


D(l) - Angle measured in the YZ plane from the Y axis toward the 
Z axis 
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F(l) - Angle measured in the XZ plane from the X axis toward the 
Z axis 

X(l) - X transfer distance 

Y(l) - Y transfer distance 

Z(l) - Z transfer distance 

NOTE: The sign of the transfer distance is taken as positive 

when the direction from the body moment reference center to the seg- 
ment centroid is in the positive axis direction 

SLOC(l) - Local reference area of the segment 

KKK(l) - This is the value of the octant in which the unit 
normal lies and ranges from 1-8. 

The following should be noted: 

1. Angles E(l), D(l), and F(l) must be between 0° and 90° • It 
should be noted that only 2 of the 3 angles are needed. If one cannot 
be determined, set that angle equal to 400. 0° and the program will 
calculate the value internally. 

2. Format - (7F10.0;I2) 

3. The octant of the unit normal may be defined as follows: 


Octant 

= 

1 

if 

unit 

normal 

lies 

in 

+X, 

+Y, 

+Z 

Octant 

= 

2 

if 

unit 

normal 

lies 

in 

-x. 

+Y, 

+Z 

Octant 

= 

3 

if 

unit 

normal 

lies 

in 

-X, 

-Y, 

+z 

Octant 

= 

4 

if 

un t 

normal 

lies 

in 

+X, 

-Y, 

+z 

Octant 

e 

5 

if 

unit 

normal 

lies 

in 

+x. 

+Y, 

-z 

Octant 

= 

6 

if 

unit 

normal 

lies 

in 

-X, 

+Y, 

-z 

Octant 

= 

7 

if 

unit 

normal 

lies 

in 

-X, 

-Y, 

-z 

Octant 

= 

8 

if 

unit 

normal 

lies 

in 

+X, 

-Y, 

-z 


Step 6 ALFHT, BETAT, FHIT Format 3F10.0 

ALFHT = a at which aerodynamics are desired 

BETAT = 3 at which aerodynamics are desired 

FHIT = 0 at which aerodynamics are desired 

It should be noted that each card counts as one toward calculating 

JA in step 4. 
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If JL equals zero, this ends the input; however, if JI has a value, 
the following steps are needed: 

Step 7 KK, Jl, J2, J*, Jb, J5, J6 Format 71^- 

KK equals 0 if only six or less segments are being changed 
equals 1 if more than six segments are being changed 
Jl - identification number of first segment being changed 
J2 - identification number of second segment being changed 
J3 - identification number of third segment being changed 
Jb - identification number of fourth segment being changed 
J5 - identification number of fifth segment being changed 
J 6 - identification number of sixth segment being changed 

In step 5* the first flat plate read in is given an internal 
identification number of 1 and each succeeding plate an identification 
number one greater until all values are read in. Therefore, to find 
the correct identification for this step, the user need only go to 
step 5 and determine what plate is being changed and its identification 
number. 

NOTE: If only one segment is being changed, J2 through J 6 are 

left blank. 


Step 8 

K(J1) D(J1) F(J1) X(J1) Y(J1) Z(J1) SLOC(jl) KKK(Jl) 


e(j6) d(j6) f(j6) x(j6) y(j6) z(j6) sloc(j6) kkk(j6) 

NOTE: If, for example, only two changes are made only 2 cards 

would be read in in this step. 



It should be noted that for all changes which have been made in 
steps 7 and 8, the program will then repeat all a, (3, and (f> variations 
that were input in step 6. 


i 


An example output is shown in table II. The first portion of the 
output consists of the title of the case and is followed by the calculated 
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direction cosines and the input quantities for each segment, the transfer 
distances, the segment area, and the angles which the normal to the seg- 
ment make. These are printed out in the order in which they were input. 
The calculated aerodynamic coefficients consisting of the force coef- 
ficients C., C„. and C %t and the moment coefficients C, , C , and C 
A 7 Y N l’ m 7 n 

are presented as C A , C^, C^ r , and C-ROLL, C-PITCH, and C-YAW. Also pre- 
sented are the aerodynamic force and moment derivatives with respect to 
sideslip (p) and angle of attack (a). 







TABLE II. - HYPERSONIC AERODYNAMIC COEFFICIENTS FOR GENERAL BODY 
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Figure 1. - Orientation of a unit normal vector of a flat plate segment 
to the principal axes and principal planes. 
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Figure 5* - Example configuration. 









